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A new class of domain wall fermions is defined that interpolates between Shamir's and Borigi's form without 
increasing the number of Dirac applications per CG iteration. This class represents a full (real) Mobius trans- 
formation of the Wilson kernel. Simulations on quenched Wilson lattices with /3 — 6.0 show that the number of 
lattice sites (Ls) in the fifth dimension can be reduced by a factor of 2 or more at fixed value of chiral symmetry 
violations measured by the residual mass (mres). 



1. INTRODUCTION 

Chiral fermions based on the Ginsparg- Wilson 
relations give a rigorous approach to exact chi- 
ral symmetry at finite lattice spacing. However 
this clear theoretical advance requires efficient al- 
gorithms to allow them to be used routinely in 
lattice gauge theory simulations. The two ma- 
jor implementations base on the overlap form of 
Neuberger ^ or the domain wall of Kaplan 
adapted by Shamir and Borigi 0] can both be 
seen as a means to satisfy the zero quark mass 
Ginsparg- Wilson relation, 

75^o.(0)-fi^o.(0)75-2i^o.(0)75Co.(0) = 275AL,, 

where the error A^^ — > is made as small as 
needed by taking Ls ~* 00, where parameter- 
izes the order of a rational polynomial approxi- 
mation. A general solution, including a non-zero 
quark mass m, is provided by the overlap opera- 
tor, 

1 + m 1 
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where e^^ [x] is the rational approximation of the 
sign function, t[x\ = x/\x\ and A^Jx] = (1 — 
e\ [x\)/A. In fact the GW relation defines an in- 
finite class of chiral fermions depending on the 
selection of the kernel D{Mr,)^ the simplest but 
certainly not the ideal example being the Wilson 
fermion operator itself. 



(1 - i^)U^,{x)5^+^,^y + (1 + i^)Ul{y)5^^y+^, 



with mass M5 = 0(— 1). On the other hand the 
standard Shamir implementation of Domain Wall 
fermions is equivalent to an overlap kernel, 
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2 + a5D^{M5) 
Here we suggest a modification of the stan- 
dard Shamir domain wall approach that replaces 
the overlap kernel, or 5-th time Hamiltonian 
75 /^(Ms), by a real Mobius transformation of the 
Wilson operator: Dw{M^) (a + hDw{-l) / {c + 
dDw{—\)) or equivalently in a more convenient 
parameterization for the domain wall formula- 
tion. 
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In addition to the Shamir parameter, = b^ — c^, 
there is a new independent scahng factor, a = 
C5 + &5, which turns out to be a significant ad- 
vantage at finite Lg even though it is irrelevant 
at Ls = 00 due scale invariance of the sign func- 
tion: e[ax] = e[x]. One way to understand the 
advantage of the Mobius transform is to notice 
that it can map any 3 real values (e.g. eigen- 
values of Dw) to arbitrary points. This allows 
one to choose M5 to separate legitimate chiral 
modes from doublers and simultaneously to ad- 
just both the smallest and largest real eigenvalue 
in D{M^). Note also Mobius fermions have a con- 
tinuous path to conventional overlap kernel as im- 
plemented by Borigi's with 65 = C5 = 05. 

2. DOMAIN WALL IMPLEMENTATION 

It is very natural to implement the real Mobius 
kernel by the 5-d Domain Wall matrix, Dowijn): 



give the direct connection to the overlap propa- 
gator, 
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where = hLOiD^{A'h) + 1 and D\ 

c^LUiDj^(M^) — 1. In this discussion we simplify 
to the case of polar decomposition {oji = 1), al- 
though Zolotarev polynomials (w^ 7^ 1) can be 
implemented as well. By standard LDU decom- 
position 5 for D^wi'm-), this leads back to an 
effective overlap operator, with 

(l+i/)i./2 + (l„iJ)-L,/2 W 

given by the polar decomposition for the "Hamil- 
tonian" H = 75 ^Mobius (-^5) defined above. The 
detailed algebraic steps will be published in a 
forthcoming article. The chiral modes at the 
boundary 
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A key observation for this action (as well as 
Borici's and Chiu's examples) is that the new ver- 
sion of "gamma 5" Hermiticity for the Mobius 
domain wall requires pulling out a factor = 
Diag[D'':}\ D^^\ ■ ■ ■ , D^_^°^] as well as reflecting, 
TZ, in the fifth dimensions: 

This in turn is reflected in the proper definition 
of chiral boundary states. 

3. DEFINITION OF RESIDUAL MASS 

The chiral breaking operator A^^ (a;) defined 
above provides the breaking term in the 4-d 
Noether theorem, 

S{TpDov{m)Tp) = m?/5(75-|-75)'(/'+2(l-m)Vi75Ai^V- 

Therefore its matrix elements are the proper mea- 
sure of the approach to exact chirality for 
00. Physically relevant matrix elements should 
be sensitive to long distance (IR) physics. In the 
Domain Wall language, there is a corresponding 
breaking term in axial Ward-Takahashi identity, 

A^J°^^'^(x) = 2TOg,A"75& + 2Q,75A°Q, 

where the Qx , Qx fields lie on any plane separat- 
ing the left and right domain walls. A convenient 
measure of the breaking term (or residual mass) 
is given by 
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in the plateau region with t away from the 
sources. An alternative definition (equivalent in 
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infinite volumes) is found by summing over all 
time slices and removing the contact term in 
Eq. |S1 This expression for m^es measures a spe- 
cific matrix element of A/,^ , 
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The positive definite kernel guarantees that zero 
residual mass implies the exact Ginsparg- Wilson 
relations and unbroken Ward-Takahashi rela- 
tions. 
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4. NUMERICAL RESULTS 

In order to have the performance of Mobius (or 
Borici) fermions on a equal footing with Shamir, 
one must be able to use even-odd preconditioning. 
However the 5-d even-odd partition is impracti- 
cal because the new operator £)_ connects even 
and odd sites so that one must invert it in each 
CG step. Instead we define the even-odd pat- 
tern only on each 4-d slice, not alternating in 5th 
axis. Now the result is that even-even and odd- 
odd matrices are independent of the gauge fields 
and can be inverted analytically at negligible cost. 
Performance tests have shown that this new 4-d 
even-odd precondition results in an improvement 
factor of roughly 3 just as the 5-d even-odd had 
done for Shamir. 

We tested the chirality and convergence of the 
Mobius domain wall operator on 20 quenched lat- 
tices at /3 — 6.0 from the NERSC gauge con- 
nection. Our results were compared with the 
standard Shamir fermions for a bare quark mass 
m = 0.06, resulting in a pion mass of roughly 
0.4 (in lattice units) for Lg ~ 16. All compar- 
isons were done for standard Shamir operator 
with as = 65 — C5 = 1.0. We tuned M5 and the 
scale a. The pion masses were adjusted so that 
Lg = S coincided with Shamir at Lg = 16. For 
Ls = 8 our tests found that = —1.5 is the op- 



Figure 1. Comparison of Lg for Mobius 

vs standard Shamir. 



timal choice, which is not surprising since it sits 
at the mid point between the two Aoki phases at 
-2.2 and -0.8. We found for Lg = 8 that the op- 
timal scale factor was a = 2.3. Furthermore we 
explored rrires at larger Lg. There appears to be 
a crossover at Lg — 16 for a ~ 4. 

More detailed performance data for a variety 
of quench and unquenched lattices is needed, but 
our preliminary observations are (1) at rrireg = 
0(0.1%) suitable for HMC simulations, Lg ^ 8 
gives essentially the same m^es as conventional 
Shamir with Lg — 16, (2) at large values of 
Lg suitable for valence quarks, rrires drops very 
rapidly if we use large scale factors and (3) fi- 
nally at fixed rrireg, the number of CG iterations 
increases only modestly for larger scale factors. 
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